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ABSTRACT 

Context. Recent revisions of the determination of the solar composition have resulted in solar models in marked dis- 
agreement with helioseismic inferences. 

Aims. The effect of the composition change on the model is largely caused by the change in the opacity. Thus we wish 
to determine an intrinsic opacity change that would compensate for the revision of the composition. 
Methods. By comparing models computed with the old and revised composition we determine the required opacity 
change. Models are computed with the opacity thus modified and used as reference in helioseismic inversions to deter- 
mine the difference between the solar and model sound speed. 

Results. An opacity increase varying from around 30 per cent near the base of the convection zone to a few percent in 
the solar core results in a sound-speed profile, with the revised composition, which is essentially indistinguishable from 
the original solar model. As a function of the logarithm of temperature this is well represented by a simple cubic fit. 
The physical realism of such a change remains debatable, however. 
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1. Introduction 

The opacity in the solar interior, and hence the solar in- 
ternal structure, depends sen sitively on the abu n dances 
of the heavy element s (e.g ., ITurck-Chieze et all . l2001at 
iTurck-Chieze fc Talonl . I2008D . Recent analyses of the so- 
lar spectrum have led to substantial revisions of the so- 
lar abundances, partic ularly of oxygen , carbon and nitro- 
gen (for a review, see lAsplundl . I2005D . Relative to previ- 
ous work these studies have the advantage of being based 
on three-dimensional hydrodynamical models of the solar 
atmosphere and taking departures from local thermody- 
namic equilibrium into account. Also, unlike earlier analy- 
ses they result in consistent abundance determinations from 
different spectral lines. As a result of the revision, the ra- 
tio Z s /X s between the present solar surface abundances 
by mass of heavy elements and hydrogen is determined 
to be 0.0165, corresponding, in calibrated solar models, to 
Z s = 0.0 125. For comparis o n, the commonly used compo- 
sition of lGrevesse fc Noelsl (|1993l ) yields Z s /X s = 0.0245, 

resulting in Z s = 0.0181. 

As pointed o ut, e .g., by Basu fe Antial 
Montalban et alj (l200l ITurck-Chieze et ail 



2001 . 



200l 



Bahcall et all (|2005l ) and iDelahave fc Pinsonneaurd T2006) 
this revision has substantial effects on solar models, 
greatly increasing the difference between their internal 
sound speed and the sound speed inferred from he- 
lioseismology. As an example, we consider Model S of 



IChristensen-D alsgaard et alj ([1996), using t he Grevesse fc 
Noels composition and OPAL opacities from llglesias et alj 
(|l992f ). and a corresponding model base d on the new abun- 
danc es and updated OPAL opacities (llglesias fc Rogers! 
I1996D . Fig. [TJ shows sound-speed differences between the 
Sun and these two models, inferred t hrough inversion o f 
the 'Best Set' of observed frequencies of lBasu "eTaLl ((1991 . 
combining data obtained with the BISON network and the 
LO WL instrument (for further details o n the inversion, 
see IChristensen-Dalsgaard fc Di Maurcl [2007). Large 
differences are also found between models based on the 
revised composition and the helioseismically inferred depth 
of the convec ti on zo ne and envelope helium abundance. 
iChaplin et alj (|2007f ) found that analysis of low-degree 
solar oscillations strongly supported the old heavy-element 
abun dance; on the ot h er han d, the results of the analysis 
by iHoudek fc Goughl (|2007[ ). with careful inclusion of 
the influence of the outer layers of the Sun, indicated a 
heavy-element abundance somewhat lower than the Model 
S value, although substantially above the value obtained 
with the revised abundances. The effects on solar models 
of the new compos ition, as teste d with h elioseismology , 
were reviewed by iBasu fc Antial (|2008l ). iGuzikl {2006) 
provided an overview of the, largely unsuccessful, attempts 
to modify the assumptions in the model computation to 
compensate for the composition change. 



J. Christensen-Dalsgaard et al.: Opacity change compensating for revised solar composition (RN) 



-o.oioh 

0.0 



0.4 0.6 

r/R 




Fig. 1. Inferred relative differences in squared sound 
speed, in the sense (Sun) — (model), fro m inversion 
of the 'Best Set' of observed frequencies of iBasu et all 
(ll997f). The open circles s h ow re sults for Model S of 
IChristensen-Dalsgaard et aTl (|T996), using the old solar 
composition, while the filled circles show results for a cor- 
responding model, assuming the revised composition. The 
horizontal bars indicate the resolution of the inversion while 
the vertical bars (hardly visible on this scale) show the 1— a 
errors in the inferences. 

By far the most important effect on solar modelling 
of the he avy-element abundance arises through the opac- 
ity. Thus [Montalban et al.l (|2004| ) found that a substantial 
opacity increase near the base of the convection zone would 
help reducing the discrepancy ca used by the revised a bun- 
dances. Similarly, it was noted bv lBahcall et al.l (|2005| ) that 
an intrinsic change in the opacity could be used to correct 
the model computation, and they estimated that an opac- 
ity increase of around 11 per cent over a relatively broad 
range in temperature would be required. Indeed, there are 
undoubtedly significant uncertainties in the very complex 
opacity calculations. In the present note, following Bahcall 
et al., we make a more detailed analysis of this nature, 
estimating the intrinsic change in the opacity required to 
obtain a model structure corresponding largely to Model S, 
but with the revised composition. 

2. Determination of the opacity change 

The goal is to determine an opacity modification such 
that the sound-speed structure of Model S can be approx- 
imately reproduced with the revised solar surface compo- 
sition. With just this constraint we can clearly only deter- 
mine a modification that depends on a single variable which 
we take to be temperature T. Thus we write the modified 
opacity as 

logK(p,T,X,Z) = log K ( As P)(p,T,X,Z) + /(logT) , (1) 

where p is density, X and Z are the abundances by mass 
of hydrogen and heavy elements, and the opac- 

ity evaluated with the revised heavy-element composition; 
log denotes logarithm to base 10. The goal is to determine 
/(logT) such that k evaluated for a structure correspond- 
ing to Model S, but with the revised heavy-element com- 
position, matches the opacity Mo del S, evaluated 
with the original iGrevesse fc Noeli (|l993f) heavy-element 
composition. 



Fig. 2. The solid curve shows the intrinsic opacity change 
AlogK = /(logT) determined from Eq. Q, on the basis 
of models with the old and the revised composition. The 
dashed curve shows the cubic fit / apP rox to A log k, given 
by Eq. ©. 



To estimate /(logT) we consider two models of 
the present Sun: one is Model S, characterized by 
{pSi Ts, Xs, Zq}, as functions of position in the model, 
and the second is Model Asp, similarly characterized by 
{/3Asp> 2as P , ^Asp, ^Asp}- We determine the difference in 
log K, at fixed T, between the two models as 

6 T log k= \og^ GN \ Ps (T),T,X s (T),Z s (T)) (2) 
- log K ( As P) (p Asp (T), T, X Asp (T),Z Asp (T)) , 

making explicit that the models are computed with k^ gn ^ 
and K^ sp ) , respectively. Using the required property of k, 
and linearizing in model differences, we obtain 



St log k 



/(logT) + 

dlogn 
d\ogX 



dlogn 
dlogpy T X 

5t log X 



St log p 



(3) 



neglecting the effect of the different dependence of Z 
on position in the two models; here, e.g., Jrlogp = 
log ps(T) — logpA s pfT 1 ) A similar expression was obtained 
by iBahcall et alj ( 20051 ). Hence, neglecting the relatively 
modest effects of the differences in X and Z at fixed T, wc 
obtain the required opacity change as 



A log k — /(log T) ~ St log k 



dlogK 
d log p 



5 T logp . (4) 



T,X 



This procedure can be iterated, to compensate for the er- 
ror in the linearization and the neglect of the composition 
effects. In practice we have found that two iterations are 
sufficient to reach a model closely matching Model S. 

3. Results 

Figure [2] shows the resulting opacity change, largely re- 
stricted to the radiative interior which evidently is the only 
region where the change is relevant. We have computed a 
full evolution sequence assuming the revised surface com- 
position and applying this change to the opacity, and, as 
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Fig. 3. Differences between Model S', computed with the 
revised composition and the opacity change A log re illus- 
trated in Fig. [2j and Model S, in the sense (Model S') 
— (Model S): Sine 2 (continuous), i51np (short dashed), 
S\np (long dashed), <51nT (dot-dashed) and SX (triple-dot- 
dashed). Here In is natural logarithm. 

Table 1. Properties of solar models (see text for a descrip- 
tion). D cz is the depth of the convective envelope, given in 
units of the solar radius Rq, and Y c is the helium abun- 
dance in the envelope. 



Model 


Dcz/Ro 




Asp 


0.2712 


0.2286 


S 


0.2885 


0.2447 


S' 


0.2882 


0.2489 


S" 


0.2875 


0.2488 



for Model S, calibrating the model to solar luminosity and 
radius as well as to the revised Z s /X s . Differences between 
the resulting Model S' and Model S are illustrated in Fig. [3] 
It is evident that the model matches Model S very closely. 
To test the effect on the comparison with the helioseismic 
inference, Fig. [4] shows the sound-speed inversion using this 
model as reference. Clearly it matches the helioseismic re- 
sults as well as does Model S. 

The simple dependence of /(logT) on logT makes it 
natural to approximate it by a low-order polynomial. In 
Fig. [U the dashed line shows the following fit: 

/approx(logT) = 0.1298-0.1856£ + 0.1064£ 2 -0.0345£ 3 , (5) 

with £ = logT — 6.3, obtained with a least-squares fit with 
uniform weight to /(logT), for logT > 6.3. The result of 
using the modification in the evolution calculation, calibrat- 
ing the model as before, and using the resulting Model S" as 
reference in a sound-speed inversion is also shown in Fig. 3] 
This is barely distinguishable from the result of using the 
original /(logT). 

Wc finally list in Table 1 the values of the depth D cz of 
the convection zone and the envelope helium abundance Y D 
for the models considered. For comparison, the helioseis- 
mically inferred value of D cz is around 0.287 and Y e has 
been determined to be around 0.25, although with some 
sensitivity to the equat ion of state used in t he solar mod- 
elling (for a review, see lBasu fc Antial . [20081. It is evident 
that Model Asp is inconsistent with the observed values, 
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Fig. 4. Inferred relative differences in squared sound speed, 
in the s ense (Sun) — (mo del), from inversion of the 'Best 
Set' of IBasu et all (|l997h . The symbols with horizontal 
and vertical bars (see caption to Fig. [IJ show results for 
Model S', computed with the revised composition and ap- 
plying the opacity correction A log k shown in Fig. [5J The 
dashed curve shows differences for Model S", similarly com- 
puted using the fit in Eq. ©. For comparison, the dotted 
curve shows the inferred difference for the original Model S. 

whereas the remaining models are essentially in accordance 
with observations. 

4. Discussion and conclusion 

The change in log k obtained in Fig. [2] corresponds to an 
opacity change of around 30 per cent near the base of the 
solar convection zone, decreasing to a very modest level in 
the core. This immediately raises the question whether such 
a change is physically realistic. Comp arisons betw e en in- 
depende nt opacity calculation s (e.g., iBadnell et all [2005) 
(see also IBasu fc Antial . I2008L for a review) indicate that 
the precision of the opacities in the relevant temperature 
range is better than 5 per cent, far smaller than the re- 
quired change. On the other hand it seems possible that, 
although highly sophisticated, the present opacity calcu- 
lation might neglect significant effects. Kurucz (personal 
communication) has noted that the neglect of a large num- 
ber of elements of low abundance could have a significant 
effect on the Rosseland mean opacity, which is sensitive to 
even rather weak absorption in spectral bands not affected 
by lines of the more abundant elements. 

It is perhaps relevant to recall the som ewhat s imilar 
situation more than two decades ago when ISimonl |l982) 
made a plea for the reexamination of opacity calculations 
in the light of problems with the modelling of certain pul- 
sating stars; he suggested that an increase in the opac- 
ity by a substantial factor could remove the discrepancies 
betwe en models and observations. Although iMagee et al.l 
(|1984| ) claimed that such an increase would be 'incom- 
patible with atomic physics' it was in fact found i n the 
OPAL opacity ca l culat ions (e.g., Ilglesias fc Rogers! . 119911 : 
iRogers fc Iglesiast Il992ft . as a result of the inclusion of the 
effects of a large number of lines. Based on this experience 
one should perhaps be wary of excluding the possibility of 
substantial opacity modifications. 

It is interesting that already the early helioseis- 
mic sound-speed inferences bv lChristensen-Dalsgaard et al.l 
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(|1985f ) indicated a need for an increase in the then current 
opacities, as subsequ ently confirmed by th e OPAL calcula- 
tions. More recently. iTripathv et al.l (|l998t ) determined the 
opacity correction to the OPAL tables required to match 
the difference between the helioseismically inferred sound 
speed and the sound speed in Model S, using the sensi- 
tivity of sound speed to opacity change s determined by 
ITripathv fe Christensen-Dalsgaardl (|l998f) ; they found that 
a relative diff erence i n opa city of order 5 per cent was 
needed. Also, I Goughl (|2004h determined the opacity dif- 
ference between Model S and a helioseismically calibrated 
solar model, finding a relative difference of order 1.5 per 
cent. 

We have assumed that the opacity correction is a func- 
tion of temperature alone. This is obviously a gross, if un- 
avoidable, simplification which should be kept in mind if 
the correction obtained here is used for other stellar-model 
calculations. Also, we emphasize that the fit given in Eq. ([5]) 
is only valid in the range [6.3, 7.2] in logT over which it was 
obtained. Even so, it might be interesting to use an opacity 
change similar to the one obtained here in stellar compu- 
tations, such as th e isoch rone analysis of M67 presented by 
IVandenBerg et all (|2007h . 

Given the difficulties arising for solar modelling from 
the revised solar abundances it is obviously crucial to 
carry out further tests of the results. Also, the prop- 
erties of the computed atmosphere models should evi- 
dently be tested a gainst other relevant observations (e.g., 
lAvres et al.l . [2006). The complexity of the hydrodynamical 
modelling of the solar atmosphere makes independent cal- 
culations highly desirable. Thu s it is encouraging that such 
calcu lations are now under way (ISteffenl.l2007tlCaffau et ail 
l2007t ). Very recently. ICaffau et al.l (|2008l ) made a determi- 
nation of the solar oxygen abundance, resulting in a value 
intermediate between the old and new compositions con- 
sidered here. We expect that the corresponding opacity 
correction required to match Model S, obtained as in the 
present analysis, wo uld be roughly ha lf the value shown in 
Fig. [2 Interestingly, lHolwegeii(|2001l ) obtained an oxygen 
abundance consiste n t with the value found by Caffau et al.; 
iTurck-Chieze et all (|2004j) showed that this did in fact re- 
sult in a sound-speed difference relative to the helioseismic 
inferences intermediate between the results for the old and 
new compositions considered here. 

Alternative independent determinations of the abun- 
dances of the relevant elements would obviously be 

very valuab le. An interesting po ssibility, proposed b y 

iGong et all (|200lD and reviewed bv lBasu fe Antial (|2008l) . 
is to constrain the heavy-element abundance from he- 
lioseismic inference of its effect on the thermodynamic 
state, and hence the sound speed, in the solar con- 
vection zone. This is an extension of the successful 
helioseismic det ermination of the solar-envelope helium 



abun d ance (e.g., [V orontso v et all 119911 : iKosovichev et all 



1992t lAntia fe Basul . 



Goughl (119841) 
the results of 



and 
such 



19941). following t h e sug gestion of 
Dappen fc Goughl (|1986h . So far 
analyses are somewhat uncertain: 



iLin fc Dappenl (120051) found sligh t indications that the 
standard ( Grevesse fc Noelsl 1 19931) heavy-element abun- 
dance wa s too high, whereas lAntia fc Basul (120061 ) and 
ILin et all (|2007[ ) concluded that the helioseismic results 
confirmed the original abundances. Further analysis, inves- 
tigations of the effects of the uncertainties in the equation 



of state, as well as better data including modes of higher 
degree, are required to obtain more definitive results. 

The relatively limited goal of the present note is to in- 
vestigate the compensating changes to the opacity required 
by the revision in the solar heavy-element composition, in 
order to match models computed with the earlier abun- 
dance determinations. This illustrates one aspect of the 
sensitivity of the solar structure, as probed by helioseis- 
mology, to the physics of the solar interior. It is evident 
that a broader goal of hclioscismology is to understand the 
full range of solar internal microphysics and dynamics re- 
quired to obtain a model in accordance with the helioseis- 
mic inferences. As discussed above, this may involve fur- 
ther adjustments of the opacity; however, it is likely that 
other processes, such as weak mixing in t he region below the 
convection zone, must be invoked ( e .g., Brun et all . Il999t 
IChristensen-Dalsgaard fc Di Maura l2007t) . Through such 
suitable adjustments to the physics used in the modelling 
it is possible to construct seis mic solar models that matc h 
the inferred sound speed (e.g.. lTurck-Chieze et all l2001bf ): 
a more interesting question is evidently the physical basis 
for these adjustments. 
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